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ABSTRACT 

-An  exact  analysis  of  diffusion  to  microdisk  and  microriag  electrodes  is  presented  for  steady-state 
conditions  and  for  the  assumptions  of  constant  concentration  over  the  surface  of  the  disk  and  conatant 
flux  over  the  surface  of  the  tint  electrodes.  Mass  transfer  to  ring  electrodes  is  progressively  enhanced  by 
reducing  their  thickness  which  also  reduces  changes  in  concentration  with  position  over  the  surface  of  the 
ring  as  well  as  the  differences  between  the  mass  transfer  coefficient  for  the  constant  concentration  and 
constant  dux  conditions.  Thin  rings  arc  therefore  recommended  for  investigations  of  electrode  kinetics. 

INTRODUCTION 

Studies  of  fast  electrochemical  kinetics  traditionally  have  been  accomplished  by 
the  use  of  transient  and  ac  perturbation  techniques,  such  as  those  pioneered  by 
Professor  Don  Smith.  Electrochemical  measurements  at  conventionally  sized  elec¬ 
trodes,  however,  invariably -contain  the  effects  of  capacitative  charging  of  the 
electrical  double  layer.  In  addition,  ohmic  losses  at  these  electrodes  can  be  consider¬ 
able  when  high  currents  flow  in  the  cell,  such  as  those  which  occur  in  potential  step 
experiments,  or  when  the  resistance  of  the  electrolyte  solution  is  high.  There  has. 
tbktaforbrbeen  a  recent  movement  toward  the  utilization  of  ultramicroelectrodes  for 
electrochemical  measurements;  the  high  mass  transport  rates,  small  currents,  and 
low  surface  areas  which  characterize  these  devices  provide  distinct  advantages  over 
measurements  at  conventionally  sized  electrodes.  High  mass  transport  rates,  due  to 
the  essentially  spherical  diffusion  Held  surrounding  these  structures,  give  rise  to 
unchanging  concentration  profiles  in  short  times^This  facilitates  the  study  of  fast 
electrochemical  kinetics  under  steady  state  conditions.  The  small  currents  flowing  in 
the  cell  reduce  ohmic  losses,  and  allow  measurements  to  be  made  in  highly  resistive 


*  Dedicated  to  the  memory  of  Don  Smith. 
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media.  The  low  surface  areas  of  these  devices  reduce  double-layer  charging  to  very 
small  proportions.  We  discuss  in  this  paper  the  behavior  of  two  commonly  used 
ultramicroelectrode  geometries. 

Most  previous  discussions  of  the  behavior  of  microdisk  (e.g.,  see  refs.  1-10)  as 
well  as  recent  discussions  of  the  behavior  of  micronng  [11,12]  electrodes  have  been 
based  on  attempts  to  solve  and  interpret  the  time-dependent  diffusion  problem  in 
the  cylindrical  coordinate  system  (or  else  using  coordinates  derived  from  this 
system): 

3c  Dd2c  D-  3c  L  _  32e 

3/  3 r1  +  r  dr  3,*J  (1) 

The  original  motivation  was  the  examination  of  the  behavior  of  comparatively  large 
electrodes  affected  by  “edge  effects."  However,  a  major  advantage  of  the  use  of 
microelectrodes  is  the  ability  to  use  these  structures  in  the  study  of  fast  electrode 
reactions  [e.g.  13,14]  (and  of  fast  reactions  in  solution  coupled  to  electrode  reactions 
[15,16])  under  steady-state  diffusion  conditions.  As  the  algebraic  or  numerical 
calculation  of  the  time-dependent  response  in  the  cylindrical  coordinate  system  is 
demanding,  it  is  therefore  always  worthwhile  to  explore  whether  there  is  a  steady- 
state  solution  for  any  given  electrode  geometry.  This  is  certainly  the  case  for 
spherical  (or  hemispherical)  electrodes  and  one  would  expect  thi«  to  hold  »l«r»  for 
microdisks  and  microrings  since  these  electrodes  are  surrounded  by  essentially 
spherical  diffusion  zones.  In  this  paper  we  therefore  explore  the  direct  solution  of 
3  2c  3  c  3  2c 

3 r2  +  +  3 z2  "  ^ 

for  these  electrode  systems. 

MICRODISKS 

The  solution  of  analogous  problems  in  heat  conduction  and  electrical  potential 
theory  is  well  known  [e.g.  17-20]  *.  As  we  shall  use  the  known  properties  of 


Lord  Kelvin  refers  10  an  experimental  determination  of  the  ratio  of  the  capacitance  of  a  sphere  and  a 
disk  as  being  1.57  (i.e.  w/2).  We  quote  Lord  Kelvin: 

“My  authority  for  this  statement  is  the  following  entry  which  I  find  written  in  pencil  on  an  old 
memorandum  book: 

PLYMOUTH.  Mond.  July  2,  1 849 

“Sir  William  Snow  Harris  has  been  showing  me  Cavendish's  unpublished  MSS  put  in  his  hands  by  Lord 
Burlington,  and  his  work  upon  them:  a  most  valuable  mine  of  results.  1  find  already  the  capacity  of  a  disc 
(circular)  was  determined  by  Cavendish  as  1/1.57  that  of  a  sphere  of  same  radius.  Now  we  have 
f  J i' 

...  •'0  <aJ-r*)1/J  a 

capacity  of  disc-— - - - -  - 

'0  (a*-r2)l/2 

It  is  much  to  be  desired  that  those  manuscripts  of  Cavendish  should  be  published  complete:  or.  at  all 
events,  that  their  safe  keeping  and  accessibility  be  secured  to  the  world.” 
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discontinuous  integrals  of  Bessel  functions  [21]  to  derive  a  solution  for  mass  transfer 
to  thin  ring  electrodes,  we  outline  briefly  the  solution  for  a  disk,  0  <  r  <  a  main¬ 
tained  at  a  constant  concentration  c‘,  with  zero  flux  over  the  surface  r>  a  i.c. 

-0  (3) 

(4) 


F-° 

dr 


-  Ac  0  <  r  <  a 
r  >  a  r  ”  0 


where  Ac  is  the  difference  in  concentration  between  infinity  and  the  surface. 
The  solution  to  eqn.  (2)  is 


c-c*-/  f(\)  exp(-Xr)i0(Xr)  dr 


(5) 


where  the  fqnction  /(X)  must  be  chosen  to  satisfy  eqns.  (3)  and  (4);  a  term 
including  in  the  general  solution  of  eqn.  (2)  must  be  excluded  as  c  remains 

finite  for  r  -» 0.  The  integral  (5)  must  be  discontinuous  at  r  —  a  while  eqn.  (4) 
requires  that 


(|f)  -  r\f(\)Jo(\r)d\-0 

\  OZ  /  *-0  •'o 


(6) 


for  r  >  a  while  for  r  <  a  the  gradient  is  arbitrary  except  that  it  cannot  be  zero. 
The  discontinuous  integrals 
r  it 

dX  I  2* 


/  sin(Xa)/0(Xr)  — 


r  <  a 

r>  a 


(7) 


associated  with 


l  1 

-oo  .  ■  r  <  a 

jT  sin(Xa)y0(Xr)  dX/  _  r2)i/2 

°  Ao  r > a 


c-c°°  -  ~(c"  ~  c‘)  j  exp(-Xr)  su>(\a)J0(\r)j£- 


satisfy  eqns.  (3)  and  (4),  respectively,  provided 
/(X)  —  (2/ir)Ac  sin(Xa) 
i.e.  we  obtain  the  solution 
2 

Ic-'-c-; 

•'o 

We  are  interested  in  the  total  flux  through  the  surface  of  the  disk 
F— 2trZ)y^|  — j  r  dr  ■»  4/)(c*  —  C>)J"  j  sin(Xo)/0(Xr)/-drdX 

-  80  /j  V  ^ 

*■  4 D(cm  -  c‘) a  f  sin(Xa)/1(Xa)-r-  -  4 Z>(c®  —  c’)a 
■'o  X 

which,  as  we  have  noted  above,  is  a  well-known  result. 


(8) 

(9) 

(10) 


(11) 
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THIN  RINGS 


It  is  most  straightforward  to  derive  expressions  for  the  concentration  distribution 
for  conditions  of  constant  flux  into  the  ring  [11,12].  The  discontinuous  integral 

oo  jO  r  >  a 

(  l/2a  r-a  (12) 

0  \  l/4a  r<a 

shows  that  a  disk  source  of  strength  Q  (mol  cm'3  s'1)  will  give  a  concentration 
distribution 

c“c"‘  +  ^dI0  exp(-Xr)/0(Xr)/,(Xfl)^  (13) 

for  the  disk  0  <  r  <  a. 

In  order  to  obtain  a  corresponding  expression  for  the  ring  a<r<b  we  suppose 
that  a  disk  sink  of  strength  -  Q  and  radius  b  is  superposed  on  the  source  (13).  This 
sink  generates  a  concentration  distribution 

c“-^of0  (I4) 

The  effect  of  the  superposed  disks  is  to  generate  a  flux  distribution  (flux  to  the 
sunace  z  -  0) 


F —  0  r  <  a  1 

F-  -q  a  <  r<  b\ 

F-  0  r>b  I 

and  the  concentration  distribution 

c”c°° +  ^rf0  exp(-*z)jo(*r)-A(*«)"^ 

exp  (-Xr)/0(Xr)7l(Xb)-y: 


(15) 


(16) 


We  are  especially  interested  in  the  concentration  distribution  over  z  -  0  as  well  as  in 
the  average  concentration  over  the  surface  of  the  ring,  c,v.  We  obtain 

Qa  r  r  ( \ / / \  _ \ Qb  /■“,/»  \ .  z , . ,  dx 


(cl)  -c~  +  ^fo  M\r)M\a)*±  -  if  f”j0(\r)jx(\b )2± 

mcm  +  %*Fi{b  -i;i;  “ }- -};i;  ~ 
Uti  Ffi.  i.  2-  fU-  l-  — 

r2  b2f  DiF'\2’  t'hb2 

QsL  F /i.  i.  2 ■  —sl\ _i£i  r/i.  l.  2-  — 

2Dr  *’  2’  r2  *i  /  22>r3/l(*’  2’  r2 


-  c"  + 


for  r  <  a 


for  a  <  r  <  b 


for  r  >  b 


(17) 


where  2F,{  }  denotes  the  relevant  hypergeometric  functions.  It  is  convenient  to 
express  eqn.  (17)  in  terms  of  a  dimensionless  parameter 

cx-o  "  (c“°  ~  C‘)D/Qb 

"I:1:  "i;1:  V1^}  f0rr<a 

f°ra<r<fc 

i:2:  i;2:  7^}  forr>6  (18) 

Similarly  we  obtain  for  the  average  concentration  over  the  surface  of  the  ring 

+  lJ Wkr>-M*'0'<l'’X 

Cv-Xr>MU>''i’f 


20a2  b 


3  irZ>(frV-fl2)  3ir£>(h2-a2)  D(b2-a2) 


2^  i;-i;  2;  ?- 


DISCUSSION 

The  expressions  derived  here  are  exact:  for  example,  eqn.  (19)  holds  for  all 
possible  values  of  b/a  and  the  known  result  for  a  disk  119] 

c„-c~  -&Qb/3wD  (20) 

is  recovered  on  setting  a  ■  0.  Furthermore,  the  expression  of  the  results  in  terms  of 
hypergeometric  functions  is  well  suited  to  computer-based  calculations.  Figure  1 
illustrates  the  relative  concentration  changes  (expressed  by  the  dimensionless 
parameter)  over  the  range  0  <  r  <  3b  at  z  -  0  for  a  disk  (b/a  -  oo)  and  for  rings 
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Fig.  1.  Plot  of  (he  dimensionless  concentration  parameter  (c*  -  c')Dfab  /as  a  function  of  the  radial 
position  for  a  disk  <f>/o-«)  and  for  rings  having  the  radius  ratios  'i'/a  —  2  and  1.2,  denotes  the 
dimensions  of  the  ring. 


having  b/a  —  2.0  and  1.2,  while  Fig.  2  illustrates  the  concentration  changes  over  me 
ring  alone.  The  figures  show  that  the  concentration  changes  are  markedly  reduced  as 
the  thickness  of  the  tings  is  decreased,  a  consequence  of  the  enhanced  mass  transfer 
to  the  ring  with  decreasing  ring  thickness.  At  this  time  it  is  certainly  feasible  to 
prepare  rings  m  10  nm  thick  on  fibers  of  10  pm  radius  giving  b/a  ■  1.001.  Figure  2 
shows  the  consequent  large  increase  of  mass  transfer  to  the  ring.  Decreasing  ring 
thickness  also  leads  to  progressively  more  uniform  concentrations  over  the  surface 
of  the  ring  which  for  b/a  —  1.001  shows  a  maximum  change  of  *  5%  between  the 
edges  and  the  center. 

Mass  transfer  coefficients  to  the  surface  of  the  ring  can  be  derived  in  a  variety  of 
ways  depending  on  the  position  at  which  the  concentration  is  defined.  A  convenient 
choice  is  the  average  concentration  [19]  and  setting  cAv  —  0  we  obtain 


km 


.2. 


P(b2-a2) 


(21) 


The  enhancement  of  mass  transfer  to  a  thin  ring  as  compared  to  that  to  a  disk  of  the 
same  radius  is  most  conveniently  expressed  by  the  ratio 


(  )ring 

i  8-K-') 

1  <*U2 

(  ^  m  )<iuk  , 

3* 

—  1 

3  it  1 

la3  ) 

1  lb  rl  1.  1  ,  “’I 

I'T’S*- ~!’2,  F/J 

where  the  suffix  Q  denotes  that  we  are  considering  conditions  of  constant  flux  over 
the  surface  of  the  electrode.  Some  values  are  listed  in  Table  1.  An  expression  for  the 


10 


b/a-ao 


0 


b 

Fig.  2.  Plot  of  the  dimensioiikia  concentration  parameter*  (cm  -  c‘)  0/ ablu  a  function  of  the  radul 
position  for  a  disk  ( b/a-ao )  and  over  the  surface  of  ring*  having  the  rfcdiui  ratio*  6/a  -  2,  1.2  and 
1.001. 


average  concentration  change  for  constant  flux  conditions  has  also  been  given  [12] 
from  which  we  derive 

(  )nn* 

(^mlduk 

■V  b>"Y'  \l-a/b)l 

The  values  calculated  from  eqq,  (23),  Table  1,  are  remarkably  close  to  the  exact 
•values  even  for  relatively  high  h/aTnshouid  be  noted,  however,  that  the  definition 
of  km  in  terms  of  the  average  concentration  implies  that  when  c^V  0  some  parts  of 
the  ring  (or  disk)  will  be  at  negative  concentrations  which  is physically  unsound. 
Alternative  definitions  based  on  the  attainment  of  c!  -  0  at  specific  positions  are 
illustrated  in  Table  1  by  the  values  for  r  —  b  and  r  —  (a  +  b)/ 2,  columns  4  and  5. 
These  values  are  again  exact  i.e.,  there  is  no  mathematical  approximation  but  the 
differences  between  columns  2,  4  and  5  illustrate  the  difficulties  of  defining  an 
appropriate  model  (see  further  below). 
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it i  s/aH  •»•«/» U 


(23) 


TABLE  1 


The  r»uo 


(< 


a  function  of  b/a  and  baaed  on 


(  ka )nng  ) 

concentration*;  comparison  wl'h  values  of  - -  ;  ~ 

(  *„)disk  /<• 


various  definitions  of  the  surface 


b/a  (*^)nng 


(*m)disk 


X0! 

approve 

Cr-, «  *»>/?  "  0 

eon.  (24) 

c*  *  const. 

cf.  ref.  22 

59°- iv25'  . 

00 

1.000 

oot  valid 

2 

1.3514 

1  3039 

1.6733 

1.284 

1.3080 

1.2747 

1.5 

1.7859 

1.7795 

2.126 

1.710 

1.7089 

1  6892 

1.25 

2.6259 

2.6229 

3.017 

2.528 

1.992 

2.4933 

2.4824 

1.2 

3.0303 

3  0279 

3  448 

2.922 

2.319 

2.8721 

2.8633 

1.125 

4.1976 

4.1962 

4.694 

4.059 

3.586 

3.9671 

3.9614 

1091 

5.3110 

5  3100 

5.833 

5.145 

4.777 

5.0112 

5.0078 

1.021 

17.396 

17.347 

18.721 

16  956 

19  372 

16.376 

16.345 

1.010 

32.887 

32.887 

40.242 

30.858 

1.005 

60.358 

60.358 

80.086 

56.573 

1  001 

254.52 

254.52 

398.84 

238.08 

The  data  in  column  4  should  correspond  for  low  (b  —  a)  to  values  derived  by  our 
previous  approximate  analysis  of  the  problem  [11];  in  that  analysis  a  uniform  ring 
sink  of  strength  Q(b-  a)  was  placed  at  r  —  (a  +  b)/2  and  the  concentration  was 
evaluated  at  r  -  b  (or  a).  This  analysis  led  to  the  simple  conclusion 

km  »*  0.4694Z)/(Z»  —  a)  (24) 

However,  the  values  derived  using  eqn.  (24)  show  that  this  approximation  deviates 
markedly  from  the  exact  expression  except  for  intermediate  values  of  b/a\  the  use 
of  eqn.  (24)  has  been  criticized  [12]. 

Table  1  also  gives  data  based  on  an  exact  calculation  of  the  capacitance  of  a  ring 
using  triple  integral  equations  [22];  the  data  in  column  7  therefore  correspond  to 
constant  concentration  over  the  surface  of  the  ring  (and  disk).  Column  8  gives  data 
based  on  an  approximate  analysis  of  this  problem  [23] 


(  ^  m  )nnf 

ir2 

(  ^  m  )duk 

1'  4(1-  £] 

\  bl 

l  U -a/bl) 

This  theory  has  usually  been  regarded  as  being  exact  [11,12,24];  Table  1  shows  that 
eqn.  (23)  is  indeed  remarkably  accurate  bearing  in  mind  the  nature  of  the  approxi¬ 
mations  made.  The  deviations  between  the  exact  values  for  the  exact  analyses  of  the 
problem  using  the  constant  flux  and  constant  concentration  approximations,  col¬ 
umns  2  and  7,  illustrate  the  difficulties  of  choosing  an  appropriate  model.  The 


constant  flux  approximation  will  apply  to  measurements  at  low  overpotentials 
(especially  for  reactions  which  are  relatively  irreversible  at  these  high  rates  of  mass 
transfer)  while  the  constant  concentration  approximation  will  hold  at  high  over- 
potentials  in  the  region  of  the  limiting  current  plateau  [11].  It  is  clear  from  the  data 
in  Figs.  1  and  2  as  well  as  Table  1  that  the  best  strategy  at  this  stage  of  development 
of  the  subject  is  to  use  micronng  electrodes  having  the  smallest  attainable  thiclcness 
since  this  will  minimize  the  differences  between  the  two  approximations. 

Although  as  we  have  shown  previously  [14],  it  is  possible  to  measure  fast  electron 
transfer  rate  information  in  the  steady  state,  it  is  clear  that  further  enhancement  in 
these  measurements  could  be  accomplished  by  using  ac  perturbation  techniques; 
this  work  will  be  the  subject  of  a  forthcoming  report. 
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